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We study the propagation of gravitational waves in a collisionless plasma with an external 
magnetic field parallel to the direction of propagation. Due to resonant interaction with the plasma 
particles the gravitational wave experiences cyclotron damping or growth, the latter case being 
possible if the distribution function for any of the particle species deviates from thermodynamical 
equilibrium. Furthermore, we examine how the damping and dispersion depends on temperature 
and on the ratio between the cyclotron- and gravitational wave frequency. The presence of the 
magnetic field leads to different dispersion relations for different polarizations, which in turn imply 
Faraday rotation of gravitational waves. 
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I. INTRODUCTION 

The propagation of weak gravitational waves in the presence of matter and electromagnetic fields has been considered 
by several authors, e.g. [Q- |Q. The back-reaction on the gravitational waves, resulting in damping and/or dispersion, 
has been studied by Refs.^J- but the gravitational effects of matter (most frequently assumed to be a fluid 
or neutral gas) and the effects of electromagnetic fields have mainly been treated separately. The interaction of 
gravitational waves with a plasma - which is the most common state of matter relevant for gravitational wave 
propagation - has been considered in an astrophysical as well as an cosmological context, see e.g. Refs. |jll|] and 
Ref. iQ respectively, and references therein. In general it turns out that matter which is in a plasma state - and 
thus exhibits electromagnetic properties - has possibilities of more efficient interaction with gravitational radiation, 
as compared to neutral matter. 

In this paper we study the propagation of weak gravitational waves in a coUisionless plasma with an external static 
and homogenous magnetic field, parallel to the direction of propagation. Naturally, the matter and fields produces a 
background curvature, but nevertheless it is meaningful to treat the background as Minkowski, provided the wavelength 
is much shorter than the background curvature (see the Appendix for a detailed discussion). It turns out that there 
is a new effect on the gravitational waves that appears due to presence of the external magnetic field - a gravitational 
analog of cyclotron damping of electromagnetic waves (cf. Ref. ||l^). In the electromagnetic case, the waves may 
interact resonantly with the gyrating motion of the particles, and the resonance occurs for particles that experience 
a wave whose Doppler shifted frequency equals the gyrofrequency. In the gravitational analog, the resonance occurs 
when the (Doppler shifted) wave frequency is twice the gyrofrequency. A mechanism for gravitational wave damping 
similar to cyclotron damping has been considered previously - namely Landau damping However, cyclotron 
damping is a potentially more important mechanism, since for this case, efficient wave-particle interaction may take 
place without the presence of ultra relativistic particles. The possibility for cyclotron damping of gravitational waves 
has been considered by Ref. and recently by Ref. 10. However, the calculation of the damping coefficient in 
Ref. do not rest on a selfconsistent kinetic theory, and it turns out that our value of the damping based on the 
Vlasov equation is smaller than theirs by several orders of magnitude. On the other hand, if we take the limit of a 
non-relativistic Maxwellian distribution function, our results essentially (see note jl^) reduces to those of Ref. |^. 

In the case where the unperturbed distribution function of the particles is not in thermodynamical equilibrium, 
gravitational wave instabilities rather than damping may occur. We give a condition on the distribution function 
for instabilities to develope and demonstrate that it may be fulfilled, for instance, by plasmas with a temperature 
anisotropy. Naturally, the gravitational wave dispersion is also modified by the presence of the magnetic field. The 
dependence of the damping and dispersion on the temperature and on the cyclotron- and gravitational wave frequencies 
is investigated. Furthermore, we confirm that the natural wave modes are circularly polarized waves also in the 
relativistic regime and that - as a consequence of the different dispersion relations for these modes (in an electron-ion 
type of plasma) - gravitational waves experiences a phenomena analogous to Faraday rotation of electromagnetic 
waves in such a medium. 
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The paper is organized as follows: In section II we present the equations governing our system, using a tetrad 
description in order to make the interpretation of our results more straightforward. Section III reviews the problem 
of single particle motion in the presence of a gravitational wave propagating parallel to an external magnetic field. 
In particular it is shown that resonant particles experiences continuous acceleration with velocities approaching the 
speed of light in vacuum. In section IV the interaction between the gravitational waves and the plasma is studied 
selfconsistently in the linearized approximation using the Einstein-Maxwell- Vlasov system of equations. The damping 
and dispersion of the gravitational wave is studied in some detail for an electron-ion plasma as well as for an electron- 
positron plasma in thermodynamical equilibrium. The case of a plasma which is not in thermodynamical equilibrium 
is studied in section V. Finally, in section VI, we summarize our results and discuss their implications. There is also 
an appendix where , the problem of separating effects into background curvature effects ("indirect" matter effect) and 
"direct" matter effects is discussed. The results in this paper, obtained using a tetrad frame formalism, are compared 
with those of a coordinate frame formalism. 



II. BASIC EQUATIONS 

We consider the interaction between weak gravitational waves and a collisionless plasma in an external magnetic 
field. Since we consider nonempty space the background space-time is necessarily curved. However, if the wavelength 
of the gravitational waves and the interaction region is small relative to the background curvature we may take the 
background to be flat and static, and the energy-momentum tensor to be the one corresponding to the perturbations 
of the electromagnetic and material fields (see the Appendix). 

For simplicity, the direction of propagation is assumed parallel to the magnetic field which we take to be static and 
homogeneous. Linearized, the Einstein field equations (EFE) take the form 

□ hab = -2k [STab - ^STrjab] (1) 

provided the gauge condition h"^^ = is fulfilled, which is equivalent to state that only tensorial perturbations are 
present. That the gauge condition is indeed satisfied will be verified in Section III below. Here □ = [c~^df — 9^] , 
hab is the small deviation from the Minkowski background metric, i.e. gab = Vab + hab, k = SttG/c^, STab is the part of 
the energy-momentum tensor containing small electromagnetic and material field perturbations associated with the 
gravitational waves and ST = ST'^. In the following it is understood that we neglect contributions of second order and 
higher in hab- In our notations a, b, c, ... = 0, 1, 2, 3 and i, j, fc, ... — 1, 2, 3 and the metric has the signature ( — h -I--I-). 

In vacuum, a linearized gravitational wave can be transformed into the transverse and traceless (TT) gauge. Then 
we have the following line-element and corresponding orthonormal frame basis 

ds2 = -c^dt^ + (1 + h+iO) dx^ + (1 - h+iO) dy^ 

+2/ix(0d.Tdy + dz2 , (2) 

eo = c~^dt , ei = (1 - \h+)dx - \h-^dy , 

62 EE (1 + \h+)dy - \h-^dx , 63 EE ^2 . 

where ^ = z — ci and /i+,ft-x <C 1. As it turns out, the gravitational waves takes this form also in the particular 
case (propagation parallel to the magnetic field) we are considering. The difference to the vacuum case will be that 
^ = z — Uphi, where Uph is the phase velocity of the gravitational wave. Fi-om now on we will refer to tetrad components 
rather than coordinate components. 

We follow the approach presented in Q for handling gravitational effects on the electromagnetic and material fields. 
Suppose an observer moves with 4- velocity This observer will measure the electric and magnetic fields Ea = Fabu'' 
and Ba = \ ^abcF'"^ , respectively, where Fab is the electromagnetic field tensor and eabc is the volume element on 
hyper-surfaces orthogonal to m". It is convenient to introduce a 3-vector notation E = (i?*) — {E^ , E'^ , E'^) etc. and 
V = e^. From now on we will assume that = c is the only nonzero component of u'^. Generally, the Maxwell 
equations contain terms coupling the electromagnetic field to the gravitational radiation field. If the gravitational 
waves propagate parallel to a magnetic field there are no gravitationally induced effects on E and B. This can be 
deduced as follows: Given the Ricci rotation coefficients for gravitational waves in the TT-gauge, the "gravitational 
source terms" in the Maxwell equations in Ref. vanishes for the given orientation of the magnetic field. 

The equation of motion for a particle of mass m and charge q in an electromagnetic and gravitational wave field is 

^p = g[E-K7m)-ipxB] -G (3) 
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where 7 = -^/l + PiP^ / {mc)^ and the four-momenta is p° = jmdx°- /dt. The gravitational force Hke term 
^^abP'^P^ /^^^j where r*^^ are the Ricci rotation coefHcients, becomes 

Gi ^{vph - Pz/im) h+pi + AxP2 
G2 = \{vph- Pz/im) -h+p2 + hxPi 



G3 



i(7m) ^ h+{pl - pI) + 2h^pip2 



(4) 
(5) 
(6) 



for weak gravitational waves propagating in the z-direction in Minkowski space, where the dot symbol represents 
derivative with respect to ^. 

In order to account for resonant wave-particle interactions we apply kinetic plasma theory, representing each particle 
species by a distribution function / governed by the Vlasov equation. In tetrad form the Vlasov equation reads 1 16 | 



Cf = Q 



where the Liouville operator is 



C = dt + {c/p'')fe, + 



and the electromagnetic force responsible for geodesic deviation is i^^M = '?(-^* + ^^'''^PjBk/irn)- In vector notation 
the Vlasov equation reads 



dtf 



pV/ 

7m 



E 



p X B 

7m 



Vp/-0 



(7) 



where Vp = (dp^,dp^,dp^). In the absence of gravitational waves, the Vlasov equation has the following spatially 
homogeneous (thermodynamical) equilibrium solution, the Synge-Jiittner distribution p^. 



fsj = 



(8) 



where uq is the spatial particle number density, ^ = mc^ /ksT, ks is the Boltzmann constant, T the temperature 
and K2{n) is a modified Bessel function of second kind. Generally, the unperturbed static solutions to Eq. (|^) 
consistent with a homogeneous and static magnetic field, say in the z -directio n, but not necessarily in thermodynamical 
equilibrium are distribution functions / = f{p±,P3), where p± = \Jp\ ^ p\. 

Since there are no induced electromagnetic fields, the perturbed energy-momentum tensor can be written 



E 



PaPh 

m7 



(9) 



where /g is the gravitational perturbation of the distribution function and the summation is over particle species 
(p.s.). The selfconsistent set of equations governing the interaction between the gravitational waves and the plasma 
are thus the two coupled equations (0) and (0). Obviously, as we have adapted a tetrad formalism, we mean the 
tetrad equivalence of Eq.(Q) | [l8[ . 



III. CYCLOTRON RESONANCE ACCELERATION 



As we will consider effects that are due to wave-particle interaction, it might be in place to first review some results 
of single particle (test-particle) motion in the presence of gravitational waves and an external magnetic field. 

Single particle motion in gravitational wave fields has been examined by many authors, see e.g. ||l|| and references 
therein. Particles moving in a monochromatic gravitational wave field experiences periodic change in its energy and 
periodic deviation from its mean direction of propagation. If the particle motion is constrained, for instance by a 
magnetic field, the change in energy and momentum may be cumulative. We refer to this as resonant acceleration. 

We focus here on charged particles in a homogenous and static magnetic field parallel to the direction of propagation 
of gravitational waves. The single particle motion in this situation have been investigated in some detail by 
treating it as a Hamiltonian dynamical problem, and also by Most noticeable is that not only can resonant 

acceleration occur, particles can even be "trapped" in such a resonant state and experience essentially unlimited linear 
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growth in kinetic energy (linear in coordinate time) and the parallel velocity will approach the velocity of light. The 
solution we present here does not - in contrast to that of Ref. - result in an exact description of the particle 
motion but does reveal the main effect (resonant acceleration of particles) and has the advantage of being straight 
forward. It also provides an intuitively clear explanation of the mechanism and the resemblance to electromagnetic 
cyclotron resonance acceleration |^^. In the electromagnetic (vacuum) case this phase-lock situation also exists, i.e., 
unlimited resonant acceleration of "trapped" particles. 

Since there are no gravitationally induced electromagnetic fields (linear in hab), Eq. (^ becomes 

= (c^e/7)px2-G = F (10) 

where luc = qB/m. The gravitational wave is assumed monochromatic and arbitrarily polarized: h+ — h+ exp[i(fcz — 
iot)] + c.c. and — exp[i(fcz — ujt)] + c.c, where c.c. stands for complex conjugate. In this section lo and k are 
assumed real. 

We assume Wc/7 ^ w/iat so that the gravitational force in Eq. (^ is small compared to the electromagnetic one. 
Thus, the particle orbits are close to the gyrating motion in the absence of a gravitational field, and we therefore make 
the ansatz pi = exp(— iti>ci/7) + c.c. and p2 — ±^p{t) exp(— IwcVt) + C-c. where the amplitude p{t) depend 

slowly on time [i.e. \dp/dt\ ^ |((jjc/7)p|] due to the gravitational influence and ± — sgn(g). Note that the cyclotron 
period is not itself a constant, but is depending on the gamma factor that also is varying slowly. Consider now the 
explicit form of the parallel driving force 



^^3 



h+p^ ± ihxP^ 
h+p Ti/ixP 



e^(kz-^t+^u^ct/l) ^C.C. (11) 



where the star symbol denotes complex conjugate. A particle with trajectory z{t) will typically experience an irregular 
oscillatory force. Unless the particle is resonant (or almost resonant) with the wave, the parallel motion will be 
random and there will be no net effect - except for the possibility of (small) diffusive acceleration ||l^ . Particles 
may, however, be resonant, i.e., have a trajectory such that there will be a non-oscillatory force resulting in a lasting 
acceleration/deceleration over several time periods of gyration. Almost resonant particles will be acted on by a force 
varying in time on a time scale (depending on the magnitude of the mismatch) slower than the gravitational wave 
period. From Eq. ( pi] ) we see that there are two possibilities for a particle to be resonant and thus acted on by a 
constant force: 

k—-uj ^ = 

dt 7 



implying 



k— oj -\ = 

dt 7 



f 12 

7 - p3/mvph 

(13) 



7-p3/mWph 

respectively, where i^ph = lojk. By z we mean the time averaged trajectory and since we will only consider time- 
averaged effects, we have put = m7dz/dt here and throughout the remainder of this section. Physically, the 
particles are resonant when they see a wave whose Doppler shifted frequency is twice the gyrofrequency Wc/7. The 
factor two - not present in the electromagnetic case - is due to the fact that the driving force is quadratic in p. Note 
that the two resonance conditions can be fulfllled simultaneously only by particles that are oppositely charged. 

The question is now whether or not the resonant particles remain resonant even though they are accelerated - a 
change in ps and 7 may potentially lead to a violation of the resonance condition ( [T^ ) or ( |T^ . Clearly the resonance 
is preserved only if 7 — ps/mvph is a constant of motion. The gravitational force, Eq. ^^)-(|^), has the property 
G ■ p = 'ymvphGs |Q. Thus it holds that 

i(^-j^]^^-J:^^o (14) 
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if and only if Vp\i = ±c, i.e. the resonance is preserved in vacuum but generally not in a medium. In many situations 
for gravitational waves, however, the vacuum relation Vph = ±c holds to a very good approximation. 

In the remainder of this section we will confine ourself to the case Wph = c, sgn{q) = — 1 and to the resonance 
condition (|l3|). We define the constant of motion a = 7 — p^/mc and observe that — mc{'-f — a) and \p\ = 
mc\/27a — 1 — a^. By considering the (time averaged) time evolution of 7 we find 

^ = w [2a-7"^(l + a2)] \h+\sm(p + \hx\ cos^ (15) 

where (p = a,rg{h+p*'^ / jm) and -0 = arg{hxP*'^ /jm). Apparently the kinetic energy is a monotonously increasing 
function for particles with suitable initial phase of the gyrating motion, i.e., for 

h = \h-^- \ cosip + \hx \ sinip > 

independently of the initial magnitude of p and p^. From now on we limit ourselves to such particles. For large times 
when 7 3> 1, Eq. ( p^ ) impHes 7 cx 2ahujt. Using the expressions for p and pz in terms of 7 and a we note that 
coordinate momenta p/7 — > as i ^ 00 whereas P3/7 scales as 

TflC 

This implies the possibility of unlimited acceleration in which p^/jm approaches c on a time scale t ^ (huj)^^. Note, 
however, that the assumption Wc/7 ^ i^hab formally restricts the predicting power of these results. However, as 
demonstrated by Ref. iQ the Hnear growth of 7 (due to parallel acceleration) remain even beyond i^c/l ^ i^hab- 

The sections to follow will concentrate on linearized gravitational wave propagation, in which case the particles 
will be assumed to deviate only slightly from the unperturbed orbits. It should be noted, however, that cyclotron 
acceleration may have interesting appHcations in the vicinity of a binary pulsars close to merging. Practically speaking 
the effective distance of acceleration close to the source will be limited by effects due to a 3-D geometry. On the other 
hand, particles close to pulsars are likely to have a relativistic background temperature, in which case the resonant 
ones may be accelerated to ultrahigh energies. 



IV. CYCLOTRON DAMPING 



As seen in the preceding section, charged particles in a homogenous static magnetic field can be accelerated and 
decelerated by gravitational waves. Thus it should not be surprising that the gravitational wave will be damped - or 
be unstable and experience growth - as it propagates through a collisionless plasma. The damping (or growth) rate 
will depend on how the particles are distributed in momentum space. In order to incorporate the damping effect due 
to this resonant wave-particle interaction mechanism, we use a kinetic description of the plasma, i.e. each plasma 
component is represented by a distribution function f[x°-,p^). The gravitational waves, /i+ = /i+ exp[i(A:z — ujt)] and 
= hx exp[\{kz — ujt)] are superimposed on the Minkowski background metric (see the Appendix), and are assumed 
to be associated with a small perturbation, fa = /cexppffcz — ujt)] , of the distribution function, i.e. f = fo + fa, 
where fo is a stationary solution to the Vlasov equation (0) in the absence of gravitational waves. This means that 
the background distribution function fo is a function of p± and p^ , where p± — \/p\ + P2- 



A. Linearized Vlasov equation 



We begin by calculating the perturbed distribution function that results from a gravitational wave propagating 
parallel to an external magnetic field. Linearizing the Vlasov Eq. (Q) in h+, hx and fa gives 



, P3 

k uj 

7m 



/g + — P X B • Vp/G - G • Vp/o 

7m 



(17) 



It is convenient to change to cylindrical coordinates in momentum space, i.e. we define pi = p± cos(f> , p2 = pi_ sine 
and ps = p|[ , so that {q/'^rn)p x B • Vp — —^~^ujcd^ and Eq. (|^ thus becomes 



7TO 



Od. 



fa — it^ [h+ cos 2(/) + hx sin 20] To 



(18) 
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where 



Eq. ( p^ has the solution 



^0 



/g 



1 



Vph Op^Jo H Op /o 

7m y 7m " 



(19) 



124, 



2 /sp|| /to — 7C1J — 20;^ 
7a;JPoe-'2'^ 



1 



2 kp\\/m — 70; + 2u;, 



- i/ix) 
[h+ + i/ix 



(20) 



The occurrence of singularities in Eq. ( po| ) indicate that there is a resonant interaction. In the case = there 
are no singularities, because without magnetic field oj > ck m our approximation, and therefore Landau damping of 
gravitational waves do not occur. For any finite value of Wc, however, the expression for fa has singularities and we 
therefore expect cyclotron damping to occur. 



B. Dispersion relation 



We now solve the tetrad equivalence ||18| of EFE (|l|) . Using (|9D and ( |20[ ) it is straight forward to confirm that the 
TT-gauge is a consistent choice in our case. The two diagonal elements both reads 



C30 27r C30 



i(/i+-i;ix)/- + ^(/i++i/ix)/+ (21) 



and the two off diagonal elements gives 



r^2 ^2 ,.21 



00 27T 00 

n3 



'-{h+-ih^)I_-^{h++ih^)I+ (22) 



where we have made use of Eq. (EQ) and 



I±-Y.^ I /pi , , dp^dp\\ (23) 

P-^- -00 

together with C = ttlokc^ /m. The two equations (^l|) and (|2^ ) combines to 

- c2fc2 - A+ - iS+] {h+ + i/ix ) = (24) 
[w^ - c^k^ -A_- iBJ\ {h+ - ihy, ) = (25) 

where 

A± = Re I± and B± = Im /± 

The natural gravitational wave modes for EFE with the given source are thus the circularly polarized modes h+ + ihx 
and /i+ — i/ix- Eqs. ( p^ and ( p^ ) contain the information about the dispersion and the damping or growth (due to 
B± 7^ 0). Due to the smallness of the gravitational coupling constant, A± and B± can typically be considered small 
in the sense that the dispersion relations Eq. (p4[) and (p5|) reads 



Lu ck + A± /2uj + iB± /2lu (26) 
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Throughout the remainder of this paper we will have this approximation in mind and we will occasionally make use 
of uj ck to simplify the expressions for I± . 

The fact that A+ ^ , unless the plasma components are of equal particle masses and as long 7^ 0, implies that 
the two gravitational wave modes have slightly different phase velocities. Thus, an incident linearly polarized wave 
- being a superposition of the two circularly polarized states - will experience a polarization shift, i.e. the direction 
of linear polarization will be rotated as it propagates through the medium. In the case of electromagnetic waves this 
is known as Faraday rotation. In principle this could be an important result as the polarization of a gravitational 
wave carries valuable information about the emitting source |p2| , e.g. the inclination of the spin axis of a quadrupole 
moment source. 

The pole contribution in the integral in ( |2^ ) is dealt with in the standard fashion, i.e. by letting the contour of 
integration pass below the pole. It should be noted that this approach disregards what happens with the distribution 
function close to the singularity, and it also misses some other aspects of the damping process, see e.g. ||2^, but it is 
the simplest way to find the main effect due to the pole, and for our purposes it suffices. 



C. Equilibrium plasma 

For a plasma in a state of thermodynamical equilibrium the unperturbed relativistic expression for the distribution 
function is the Synge-Jiittner distribution /o — fsj, defined in Eq. (||), and for this choice J-"o reduces to 

2 

■^0 ^ "■^w (27) 

As no instabilities can develope the gravitational waves will exhibit damping. In this section we examine how the 
damping and dispersion depends on the ratio ujc/uj and on the temperature of the plasma. First we examine the 
nonrelativistic regime for which analytic results can be obtained. In the relativistic regime we then present results 
for an electron-positron plasma and an electron-ion plasma obtained by numerical integration of Eq. ( p3| ) . For this 
purpose it is practical to introduce normalized momenta pj^ = p±/mc, = p\\/mc and a dimensionless frequency 
ratio fl = u>c/w. 

1. Nonrelativistic temperature 

In the regime of nonrelativistic particle velocities Eq. (^) together with (^) becomes 

327r / /i \3/2 7 , 



p.s. 



after performing the pa_ integration, where we have introduced t = 1/^/TrmnoG, which is the characteristic time 
for gravitational contraction of a g as w ith density no and particle mass to. In the case of a cosmological plasma 
T coincides, apart from a factor a/3/8, with the Hubble time. From Eq. ( p^ ) we obtain to lowest order in the 
temperature 

v-^ 16 1 , , 

p.s. ^ 

_y^8V^^_^(1^2^)./2 

Note that for Eq. ( p9| ) to apply we must assume ± 2fl)^^ ^ 1, i.e. for the wave frequency close enough to 

twice the gyrofrequency, higher order thermal effects is always important for the wave dispersion, whereas Eq.(|30|) 
- which follows from the residue theorem - holds for all values of fl. Still, it is clear that the dispersion can be 
enhanced by the magnetic field due to the existence of the resonance. The exponential decrease of the damping with 
1 =F 2fJ implies that significant damping only occurs in a limited region in frequency space close to 1 =F 2ri = for 
a low temperature plasma. In the limit of zero temperature the size of this region tends to zero. The magnitude of 
both A± and B± are monotonically increasing with temperature and vanishes at zero temperature. The results in 
this subsection essentially agrees with Ref. ||^, see note In the case of no magnetic field, $7 = 0, the dispersion 
relations (HJ) and ( p5| ) both reduces to — c^k'^ — A = 0, where A = J^p s ^Q/T^f^- This expression is in agreement 
lEl with Refs. " 
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2. Electron-positron plasma 



Denote Tcp = l/VTrmcngG (e and p stands for electron and positron, respectively). In this case only the sign of the 
charge differ for the two particle species and we may write Eq. (j2^), applying (p7|), as 



/±(f^) = -r,p2V-iVT / I 1 rfPi^PlI (31) 



where = Wc/^^ and a)c = \q\B/m. Note that /±(ri) = = I±{—^)- The normalized function T±(Cl) = 

-f±M'^cp is numerically calculated for four different temperatures, namely ^ = 100 ("nonrelativistic") 10, 1 and 0.1 
("ultrarelativistic") . The results are displayed in Figure 1 (note that T± is normalized against the temperature). In 
the nonrelativistic case (see Figure 1. (a)) we have a finite region, —1/2 < < 1/2 , where Re/± is positive. Re/± 
changes sign at J7 « ±1/2 and approaches zero like 1/ (t^) in the limit ±oo. Im/± is negative definite (and thus 
there is indeed damping of the gravitational wave for any finite fl) and has a gaussian shape around the resonances 
n « ±1/2. 

This is also the characteristic behavior for Re/± and Im/± at higher temperatures. The relativistic effects on the 
dispersion and the damping are the following: i) The resonance peaks of Im/± (which occur at fl — ±1/2 at zero 
temperature) are shifted to higher values of for higher temperatures. Also Re/± experiences a similar shift. In the 
ultrarelativistic case (see Figure 1. (d)), Re/± changes sign at 17 ~ ±25 and Im/± is centered about fl ~ ±20. ii) The 
magnitude of both Re I± and Im I± increases roughly linearly with temperature (note that the curves are normalized 
against {fJ'T^p)~^ in Figure 1.), which can be compared with the temperature dependence in the nonrelativistic regime 



that is given by Eq. (|2S|) and Eq. (|30|). Hi) The region of damping broadens. In the cold Hmit ImJ± takes the 
form of two separated gaussian functions, the width tending to zero with diminishing temperature. These regions are 
widened and the gaussian shape is deformed with increasing temperature. In the ultrarelativistic case Im/± decays 
exponentially as 17 ^ ±cxd but approaches zero more abruptly as 17 — + 0. 



3. Electron-ion plasma 



Denote r = Tei = l/y/TrrriinoG (e and i stands for electron and ion, respectively). For an electron-ion plasma there is 
an asymmetry between the particle species due to the small mass ratio e = rric/mi, giving different order of magnitudes 
for the two cyclotron frequencies uucc and Ud- Thus the resonances will occur for very different gravitational wave 
frequencies. Given Eq. (p7|), we have the following form of Eq. (ESh 



^M_IJ ^i + pi + p|p||-,yi + pi + p|±2u 



-rT^-^ / / , ' ^dp^dp^^ (32) 



where 



Mo/i = , rr. and n^/ 



C I 



BJ-c/i ^ "l-c/i^ 

Provided that Tq = Ti, /ij and differs typically by three orders in magnitude. This implies that the ions can be 
considered nonrelativistic even in the regime of ultrarelativistic electrons (/Zo ^ 0.1). It is convenient to consider the 
two frequency domains 17e ~ 1 and fii ~ 1 separately. Noting that 17i = el7e and in = iie, in the region lie ~ 1 we 
have 



/±(17e)«-T-2-^ / / ^J^= ' _dp^dp||+4^ '-^ (33) 

K^iiic) J J /l+p^+p2p /i+p^+p2±217, " r^p^^cl±2en, 
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The normalized function X+(f2o) = I+l^^'^p is displayed in Figure 2 for = 100 ("nonrelativistic") 10, 1 and 0.1 
("ultrarelativistic") . By symmetry, /-(Sic) is the mirror image of /+(17c), i-e. I-{^c) — I+{—^o) and therefore this 
curve is not presented. Except for the lack of symmetry about fio = the result is similar to that of the electron- 
positron plasma. In the given temperature and frequency domains the ion contribution to the normalized function 
/±(f2e) is just the approximately constant value 16. 

The effect of resonant ions becomes important in the frequency domain ~ 1. In Fig. 3 we show 2+{Vti) = I+^it^^ 
for jjLe — 0.1. The corresponding figures are qualitatively similar in the entire temperature domain 0.1 — 100, however, 
and therefore only one of them is shown. The small region near lli = (fii '--^ e to be specific) that contains the 
electron contribution has been left out. 



4- The group velocity 



The fact that A± is at some points negative implies, together with (|2q), that for some wavelengths and frequencies 
the group velocity of the waves exceeds the speed of light. Superluminal group velocities for gravitational waves 
have been found before, see Ref. |25| and references therein. In most cases, but not in all, it has been an effect of 
the background curvature. It should be noted that several results from the literature are in contradiction with each 
other. Naturally, in our case the superluminal group velocity is a direct effect of the medium. For the case of an 
electron-positron plasma the group velocity corresponding to Eq. (p^) is 



duj 
Ik 



1 

2^ 



A, 



''-A' 



(34) 



where the prime denotes derivative with respect to Cuc/^^- For ujc = 0, i.e. in the case of no magnetic field, the medium 
is just a collisionless gas of charged particles and the group velocity is smaller than the velocity of light, in agreement 
with the results of previous authors, e.g. From Figure 1 it is clear that there are regions where Vg > c is realized, 
for instance about the point Cjc/uj ~ 0.6 (where A± has a local minima) in Figure 1(a). Similarly, also an electron-ion 
plasma allows superluminal group velocities. Group velocities that exceeds the speed of light is not necessarily at 
odds with causality - an issue explored for instance in Ref. p5||. The group velocity can simply not be interpreted as 
the (gravitational wave) signal velocity in this situation. 



V. NONEQUILIBRIUM PLASMA 



In the case of thermodynamical nonequihbrium, the system has free energy that may feed a gravitational wave 
instability. This occurs whenever the imaginary part of I±, defined by Eq. (gSj), is somewhere positive. Applying the 
residue theorem, Eq. ( |2^ ) gives 



oo 





\ ^" 1 


fpl 


Aj-k/m 



dpi 



(35) 



where pr = mc (71. ± 2uJc/uj) is the resonant parallel momenta, 

2u!r 



7r 



Lu — ck 



(o^^-c^ 
4^2 I 



p1 



is the gamma factor evaluated at the resonant momenta and 

UJ Pr 



A,. = 1 



ck jj Tnc 



1/2^ 



Instabilities occurs for distribution functions and frequencies such that the condition B± > is satisfied. For simplicity 
we focus on nonrelativistic temperatures from now on, in which case Eq. (pa) reduces to 
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^± = E /^^2 2 / pI [4Pc/o + Pi5p„ /o] dp± (36) 

p.s. (, 

where Pc = ±2mccjc/w- In order to show that there indeed exist instabihties we consider the following example of a 
temperature-anisotropic "drifting Maxwellian" distribution function 

= 3,/° , e-[iv,~v.?/vl.,+vl/pl.A (37) 



where pth|| = y^^fceTji , ptii± = \/2mkBT± and pd is the drift momenta, for which B± becomes 

B^^y ae] e-™^^^«VpL„ (38) 

where a = T±/T\\ — 1 and 8 = 1 ± 2ujcl io — ■p^l'tnc. It is easily seen that B± > for certain values of a and u>c/uj. Note 
that if a = then B± < for all values of lOc/uj and hence there can be no temperature- isoiropic beam instability 
|p6| . On the other hand, putting pd = 0, it is clear that a temperature-anisotropic distribution function without a 
drift can be the source of an instability. 

The gravitational waves produced by a homogeneous plasma, due to the above cyclotron resonance instability will 
have a frequency of the order of the cyclotron frequency, at least in the nonrelativistic temperature limit considered 
here. In principle this opens up the possibility of emission of high frequency gravitational waves through large 
magnetic fields (of the order of 10~^ T or larger), i.e. frequencies well above the frequency range expected from 
"conventional" gravitational wave sources, such as compact binaries, neutron-star normal modes and gravitationally 
collapsing objects, reaching up to 10 kHz. Still the radiation considered here is generated by a collective process where 
the amount of matter interacting coherently can be as large as in other astrophysical examples. 

It is not so easy to find astrophysical applications of the cyclotron resonance instability, however: Firstly there 
must be a magnetized plasma cloud, with a nonequilibrium distribution function that fulfills the condition B± > 0. 
Such a cloud could in principle be generated, for example if there is a magnetic field geometry that allows for a loss- 
cone distribution to evolve. However, there is an obvious risk that there will be purely electromagnetic instabilities, 
which typically have much higher growth rates than the gravitational one, that will dominate the picture. Secondly, 
for significant gravitational generation, the plasma cloud must be very much denser than the average density of the 
universe, otherwise the growth of the amplitude will take place slowly even compared to a cosmological time-scale. 
In a dense plasma cloud, on the other hand, gravitation must be balanced by pressure gradients in order not to 
selfcontract before significant radiation generation due to the cyclotron resonance instability occurs. This suggests 
that maybe plasma inhomogeneities should be included in our treatment, although it seems likely that the effects 
of inhomogeneity may be neglected as long as the gradient scale lengths is much longer than the wavelength of the 
gravitational radiation. 



VI. SUMMARY AND DISCUSSION 



We have considered linearized gravitational waves in the short wavelength approximation, propagating in a plasma 
parallel to an external magnetic field. In vacuum, there is the possibility of cyclotron acceleration of charged particles 
up to velocities arbitrarily close to the speed of light. Taking the collective effects of particle distributions on the 
gravitational waves into account, it follows that the (gravitational) normal modes for the system are circularly polarized 
gravitational waves, and we derive the corresponding dispersion relations, which coincide with that of Ref. js) in the 
limit of a low-temperature Maxwellian plasma. In the case of an equilibrium plasma the waves are shown to be damped 
due to resonant interaction with the plasma particles and the dispersion is modified and enhanced as compared to 
the case of no magnetic field. In the case of thermodynamical nonequilibrium, there is the possibility of gravitational 
wave instabilities. To show in a concrete way that this can be reaHzed, we demonstrate that there are temperature 
anisotropic distribution functions that are unstable. Furthermore, we have examined how the damping and dispersion 
in an electron-positron type of plasma and an electron- ion plasma, respectively, depends on the ratio Wc/w and on T 
in the regime of relativistic temperature. The strongest effects occurs when u and ujc are comparable and the effect 
increases with temperature and density. 

The question is whether cyclotron damping can be observed, if we assume that gravitational wave astronomy ||2|] 
develops successfully. The calculations made in Ref. ||lj| end up with an estimate of 10% damping during a propagation 



10 



distance of the order of 30 kpc, which suggests that there is at least some chance of gravitational cyclotron damping 
to be observed. Combining Eq. with the definition I± = I±ht'^ (recall that X± is the normaHzed value of I± 
displayed in Figurel, 2 and 3) we note that our damping rate is 

Similarly, the time-scale for gravitational wave dispersion is /xwr^/ ReX±. For typical values for (equilibrium) plasma 
in interstellar space, it is clear that the damping rate predicted by Eq. ( ^9| ) is several orders of magnitude smaller 
than the estimation made in Ref. [Q. Presumably the discrepancy is due to the fact that in Ref. only the effect 
of acceleration of particles (corresponding to energy being transported from the gravitational wave) is considered. 
Generally when the damping due to wave-particle interaction is small, it should be noted that as wtprop (where tprop 
is the time of wave- particle interaction) grows, the energy loss of the decelerated particles becomes very close to 
the energy gained by the accelerated particles, and thus the omission of the contribution from decelerated particles 
leads to large errors in the damping coefficient. If the singularities is treated properly, the effect of decelerating as 
well as accelerating particles on the wave is included automatically in a kinetic description. In general for parallel 
propagation we find that cyclotron damping as well as dispersion of gravitational waves through interstellar space is 
neghgible, in the sense that there is essentially no hope of detecting it with gravitational wave detectors with realistic 
sensitivity. 

In general gravitational waves propagate in an angle to the magnetic field. It is well-known (see e.g. ||l||) that this 
leads to generation of electromagnetic fields by the gravitational wave, and a dispersion relation governing propagation 
perpendicular to a magnetic field in a plasma has been derived by Ref. |^. The results shows that the gravitational 
wave is most affected by the matter when the frequency matches one of the natural frequencies of the system, like 
the cyclotron frequency or the plasma frequency. It should be noted that in an inhomogeneous medium with slowly 
varying background parameters we will typically reach a point where the generated electromagnetic fields fulfills the 
dispersion relation of some natural plasma mode. In that case Hnear mode conversion, which is a resonant process 
involving all particles, may take place. However, such a problem remains an issue for future research. 



VII. APPENDIX 

At a first sight our results seem to disagree with those of Ref. |||]. However, we will show below that our results 
essentially are in agreement ||5|, although the comparison is non-trivial. In general the effect of a curved background 
space-time cannot be neglected. However, for linearized gravitational waves with short wavelength compared to the 
background curvature, the back-reaction on the gravitational wave can be separated into two effects, those respective 
contribution can be added to the fiat vacuum dispersion relation o;^ — fc^c^ = 0, see Ref. 1 p. 427. The first effect (that 
is considered by us) is a direct consequence of matter and fields (where the background curvature can be neglected) , 
and the second effect is an indirect consequence due to the background curvature produced by the matter and fields . 
How to find the curved background contribution to the gravitational wave dispersion relation for a Robertson- Walker 
or a Schwarzschild background metric, see Ref. |Q and Ref. |0| respectively. Here a warning is strongly motivated, 
however. In general the separation of contributions to the dispersion relation into the "direct" and "indirect" effect 
is not completely unique, but to some extent depends on the formalism used. Thus great care must be taken when 
adding a "direct" and an "indirect" contribution to the dispersion relation derived by different authors, to see that the 
formalisms that have produced the different expressions are compatible. 

To illustrate the above matters we compare our calculations using a tetrad frame formaHsm with those of a coordinate 
frame formalism. In the absence of gravitational wave perturbations the Einstein field equation reads 

In this section we use Greek indices /i, ly, ... — 0, 1, 2, 3 (or t, x, y, z) for coordinate components and reserve a, b, ... — 
0, 1,2,3 for tetrad components. Adding a small perturbation, so that R^i, = r!"^} + 5R^v, T^i, = T^J + 5T^i, and 
giiv = glS + hfii, , gives the linearized equation 

= K(5[T^, - y^,T] (41) 

Focusing on the direct effect of matter and assuming the short wavelength regime we use gjlS = , and after suitable 
gauge transformations we obtain: 

□ h,,, = ~2k [ST^, - ^STt]^,] (42) 
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It should be noted that in principle the left hand side of Eq. ( [4 21 ) should contain "cross terms" proportional to 
the product of /i^j/ and components of the background part of the Riemann tensor (see Ref. corresponding to 
the "indirect effects of matter" but these terms are omitted, not because they are small, but because in the short 
wavelength limit their contribution can be calculated separately and added afterwards. For the same reason also a term 
proportional to T^^^h^^, = —R'^'^^h^^, have been omitted from the right hand side. Eq. ( ^2| ) also applies to the tetrad 
description if it is understood that hab denotes what SRab reduces to in this approximation, i.e. hu = —h22 = h+ and 
hi2 = h2i = hx in the TT-gauge. The short wavelength approximation and TT-gauge is now assumed throughout 
the remainder of this section if nothing else is said. 

The total energy-momentum tensor due to electromagnetic fields and matter is 



F'^ F 



J rn-f 



211/2 



where is the magnetic permeability, 7 = pt/mc = [l + P^/ (mc)^ — h^^'''pfj_p^/{mc)^] and pt is the zero component 
of the four- momenta satisfying = —rri^c^. Note that the definition 7 = pt/mc differs from the tetrad formalism. 
The background magnetic field is taken to be F12 — F^"^ = B (being identical io that in the tetrad description). Hence 
the linearized energy-momentum tensor reads 



1 . 



E 

p.s. 



77170 



-E 

p.s. 
I 1/2 



(43) 



Note that only the last term on the right hand 
since we have gab = Tjab in that case. In some 



where ^7 = —^h^^'^p^p^/ (7omc)^ and 70 = [l-f p^/(mc)^ 
side appears in the tetrad equation corresponding to Eq. ( |43[) 
sense the perturbation of the electromagnetic part of the energy momentum tensor is somewhat artificial, since the 
electromagnetic field itself is not perturbed, although it is clear from a technical point of view that all terms in Eq. 
( ^3| ) must be included. 

In the coordinate description the Vlasov equation reads 



8J + ^ 



V/ + [F - G] • Vp/ = 



where — qp^F^^/jm, = T^^^jP^^p" /^m, F^^, are the Christoffel symbols corresponding to g^jy — -q^u + and 
V ~ {dx,dy,dz)- The Hnearized Vlasov equation for the gravitationally perturbed part of the distribution function, 
fa, thus becomes 



dt 



TO70 



• V Fo • Vp 



fa + [SF - G] • Vp/o = 



where 



(44) 

(45) 
(46) 

and JFq was defined in Eq. (|19|). The term 6F in Eq. ( |4^ ) has its origin from the lowering of an index on F'^ and 
therefore does not occur in the tetrad description [Q. Also the term G • Vp/o differs from the tetrad description 
because Gi = ^G^, G2 = ^Gy, and G3 = G^. 

Solving Eq. (p3) and applying the solution to Eq. (42) gives the dispersion relations 



SF ■ Vp/o = 2ijJc7o ^ [^+ sin 20 — hx cos 20] To 
G • Vp/o = i (2ti> — fcp|| (7om)^^) [^,4. cos 20 — hx sin 20] J^o 



where A± = Re /± + ^o , B± 



c^k^ - A+- \B+\ {h+ + i/ix) = 
[uj^ - c^P - A_ - {h+ - ihx) = 
Im/± and 



(47) 
(48) 



E 



p1 



2u} - A;p||/(7om) ± 2^^/70 
kp^^/m - 70W =F 



Todp±dpi\ 



(49) 



Ao = KcVo'5' + ^E 



p.s. 



-00 



homy 



(50) 
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A couple of things should be noted. Firstly, the coordinate frame description gives a term contained in proportional 
to that remains even in the absence of particles. This term has no correspondence in the tetrad formaHsm, and 
the reason is that there are no perturbations of the electromagnetic part of the energy momentum tensor when we 
study the projections on the tetrad basis vectors, in contrast to what happens when we look at the coordinate basis 
components (cf. Eqs. (^) and (^3|)). However, the division into the perturbed and unperturbed tensors are the 
same independent of formalism, and thus the "direct term" proportional to B"^ in the coordinate frame formaHsm 
is compensated by a term associated with the background curvature produced by the unperturbed magnetic field 
in the tetrad formalism. On the other hand the division into a direct effect of matter and an indirect (background 
curvature) effect is not completely artificial: The direct effect is determined entirely by the local matter content, 
whereas the background curvature effect is determined by the matter content both locally and globally. Furthermore, 
the background curvature effect will necessarily give contributions to ( p7| ) and ( p8| ) that are real and independent of 
the gravitational wave frequency (cf. the discussion after Eq. (E3)). Noting that 



2uj - A:p||/(7om) ± 2ujc/lo _ uj 

kp\i/m — jQUj ^ 2ujc kp\\/m — jquj =f 2uJc 



■ 70 (51) 



it is clear that the correction to the vacuum dispersion relations in the coordinate and tetrad frame formalism deviates 
only by a frequency independent real constant, and in particular the cyclotron damping agrees perfe ctly within the 
short wavelength approximation scheme. Adding the contribution to I± from the second term in Eq. ( pi] ) to Ao then 
all differences between the tetrad and coordinate frame formalisms are collected in this term: 



KC^^Q ^B^ + TTK 

p.s. 



-oo 



r/0 - 7oPj_-^oc /m 



2(7om)3- 



dp\\dpA_ (52) 



Relaxing the assumption of Minkowskian background, i.e. adding "cross terms" (cf. Ref. 0]) to Eq. ( ^ (which of 
course requires a different tetrad basis in that formalism) , we could in principle confirm the agreement of the total 
dispersions relations including both direct and indirect effects. This would be a tedious task, however, since we must 
then solve for the background configuration including the anisotropic magnetic field contribution. Since we have 
shown that the separation into background curvature and "direct" matter effects is not unique, one can question the 
relevance of the figures showing the real (dispersive) part of I± since only the direct matter effect is contained (recall 
that the damping contribution is unique, however, and due to the direct effect only). But it turns out that there 
are three reasons that make also the real plots of I± relevant: Firstly, the main contribution to the real value of 
I± occurs in the frequency regime when cyclotron resonance effects plays a role, which always is a direct effect of 
matter independent of formalism. Secondly, for many cases the dominant contribution to the background curvature 
comes from the rest mass (Too) part of the energy momentum tensor, and if that is the case, the difference displayed 
in (IH^ ) is negligible as compared to other background curvature contributions [ p7| . For example we could add the 
background curvature contributions for a Robertson- Walker or a Schwarzschild background metric, using the results 
of Ref. ||l| and Ref. |@| respectively, since the difference in formalisms used in most cases would give only a negligible 
contribution in comparison with that included. Thirdly, since the background curvature gives only a real constant 
contribution to the dispersion relations, our plots always reveal the frequency dependence of the dispersion. 
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FIG. 1. The real (solid line) and imaginary (dotted line) part of T± = I±^T^^ for an electron-positron plasma at four 
different temperatures: (a) /x = 100, (b) /x = 10 ,(c) /x = 1 and (d) jJL — Q.l. 

FIG. 2. The real (solid line) and imaginary (dotted line) part of = Ij^jjLT^^ for an electron-ion plasma at four different 
temperatures: (a) /i = 100, (b) /i = 10 ,(c) /x = 1 and (d) pi = 0.1. 

FIG. 3. The real (solid line) and imaginary (dotted line) part of X+ = I^j-iT^^ for an electron-ion plasma for jig = 0.1. The 
small region near fij = containing the electron contribution has been left out. 
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